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Three  studies  related  to  the  testing  of  large  scale  bubble  memories  were 
undertaken.  In  the  first,  a  technique  was  developed  which  optically  scanned 
a  bubble  memory  overlay  pattern  to  identify  and  locate  faults  in  the  pattern. 
Some  preliminary  experiments  indicate  the  feasibility  of  this  approach.  The 
second  study  looked  into  the  microscopic  positioning  of  the  bubble  memory  chip. 
It  was  found  that  the  positioning  mechanism  can  be  consistent  with  accuracies 
of  O.^us.  Thus  1  tun  element  overlay  patterns  can  be  scanned  as  indicated 
above.  Finally  an  integrated  optics  approach  was  examined  for  application  to 
dynamic  testing.  Although  much  work  remains  to  be  done,  it  appears  that  the 
approach  described  here  is  possible. 
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SECTION  I 


Materials  Research 


A.  M-H  Looper 
1.  Introduction 

The  characterization  of  the  static  parameters  of  a  magnetic 
bubble  material  traditionally  consists  of  three  parameters.  The  fill, 
thickness  h,  the  stripewidth  w,  at  zero  applied  field  and  the  bubble 
collapse  field  Hbc  are  measured^1}  All  pertinent  static  parameters 
can  be  calculated  from  these  three  terms.  These  include  the  saturation 
magnetization,  the  domain  wall  energy  density,  the  material  charac¬ 
teristic  length  and  others.  Ordinarily,  the  stripewidth  and  collapse 
field  measurements  are  done  optically  with  a  microscope.  As  bubble 
sizes  decrease,  the  difficulties  in  these  measurements  increase. 

A  determination  of  the  hysteresis  loop  of  a  wafer  provides  an 
alternative  to  optical  measurements.  Measurement  of  the  initial 
susceptibility,  x0  =  (dM/dH)m_Q  and  stripe  collapse  field,  H$c  from 
the  hysteresis  loop  are  substituted  for  w  and  Hbc> 

From  the  stripe  domain  stability  conditions  developed  by  Kooy 
and  Enzf2)  the  relationship  between  the  film  magnetization  and  the 
applied  magnetic  field  is  given  as: 


4ttM-H  2/yj 
4nMs  1  +  /y 


sin(irn(l  + 


[1  -  exp(-2nno)]  =  0. 


(1) 


M  is  the  net  magnetization  of  the  film,  Mg  is  the  saturation  magnet¬ 
ization,  H  is  the  applied  field,  n  is  the  film  thickness, 
v  =  1  +  2*Ms/Ku>  a  =  h/u/(d1  +  d2),  di  is  the  stripe  width  for  domains 
with  ft  is  the  direction  of  ft  and  d2  is  for  domains  with  ft  in  the  op¬ 
posite  direction.  The  equation  assumes  a  film  of  infinite  area  with 
straight  domains. 


1 


The  domain  wail  energy  can  be  related  to  the  magnetization  by  the 
equation: 


0 

T4Ms)=h7u 


4-n2a2 


_ 00 

Z  -Si 


sin2[~  (i  .  Jj-)m  . 


n=l 


(1  +  2TTna)exo(-2t!nci)  ] 


(2) 


In  this  equation,  c>w  is  the  energy  of  the  domain  wall  per  unit  area. 

If  the  uniaxial  anisotropy  field  is  large  compared  to  4wM  and 

3  s 

u  -  1,  the  above  equations  reduce  to  : 


^r1  *  ^  Z  -b  sin"”'<1  *  t)"1  -  ' 0 

s  s 

n=l 


(3) 


and 


Z  si"2((Jr)  a 

n=l 


.  exp(^li)(1  ,  2^),  (4) 


where  d  d;  +  d;  and  x  is  the  material  length  (h  o>w/4*M2). 

From  equation  3  the  initial  susceptibility  is  calculated. 


4n,  =  l1  +  ^  S  n  (*1,ni1  '  exP^)H_1  (5) 

nsl 

where  w  h  d/2.  The  relation  between  x  and  w/h  can  be  graphed  or  tab- 

0  13) 

ulated  and  if  \0  is  measured,  w/h  can  be  determined  •  Similarly,  using 
equation  4,  x0  and  x/h  can  be  related.  Using  Thieles  conditions  for 


2 


(4} 

bubble  stability  '.graphs  of  4nMs  against  stripe  collapse  field  can  be 
drawn.  If  an  independent  measurement  of  h  is  made,  w  and  4irMs  are 
determined. 

2.  Principle  of  Operation 

The  principle  of  operation  of  the  hysteresis  loop  device 
("M-H  Looper  )  is  to  place  a  sample  in  an  alternating  magnetic  drive 
field  (Hd  =  Hq  sin  wt)  and  by  using  a  small  pick  up  coil  placed  in 
close  proximity  to  the  sample,  measure  the  sample's  magnetization  as 
a  function  of  applied  field.  Figure  lisa  block  diagram  of  the  sys¬ 
tem  including  sample  positioning  and  pick  up  coil  placement. 

Since  the  change  in  the  magnetic  flux  linkage  through  the  sample 
pick  up  coil  will  be  dominated  by  the  changes  in  Hd  rather  than  the 
change  in  the  sample  magnetization,  a  second  "compensation"  pick  up 
coil  is  required.  This  second  coil  is  placed  in  a  location  such  that 
the  changes  in  can  be  monitored  while  the  change  in  the  sample 
magnetization  will  have  a  negligible  effect  on  the  induced  voltages 
in  this  coil . 

The  signal  from  the  compensation  coil,  Vcomp(t),  is  amplified 

and  subtracted  from  the  amplified  sample  coil  signal,  vsamp(t),  such 

that  V  .  =  V_  -  V  „  =  0  when  no  sample  is  present.  Thus,  when 

sig  samp  comp  r 

a  sample  is  present,  V$1-  ( t)  is  directly  proportional  to  the  time 

derivative  of  the  sample  magnetization,  i.e.  V  .  01  dM(t)/dt.  Signal 

b  i  y 

processing  and  hysteresis  loop  display  are  fully  automated  using  a 
desk  top  calculator  system.  Figures  2  and  3  are  photographs  of  the 
system.  The  desk  top  calculator  and  plotter  are  not  shown. 

a.  Drive  System 

For  operation  of  the  looper,  rapidly  varying  (for  resonable 
signal  strength)  homogeneous  fields  of  sufficient  intensity  to 
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Figure  1.  Block  diagram  of  M-H  looper. 


Figure  2.  The  M-H  looper.  Shown  are  the  frequency  generator,  the  bipolar  power  supply,  the  sense 
resistor  box,  the  capacitor  bank  and  the  drive  coil  assembly. 


Figure  2.  The  V-H  loop er.  Shown  are  the  electronics  and  co~persa tier  bridge  box,  the  oscilloscope 

and  the  calculator  controlled  .-oltreter.  7  re  sa~cle  <.a'er  'tool  ard  sense  resistor  sicnals 


saturate  bubble  materials  are  desired.  These  fields  are  produced  with 
a  set  of  Helmholtz  coils  driven  by  a  Kepco  15-20  M  bipolar  operation¬ 
al  power  supply.  The  coils  have  300  turns  of  #14  gauge  copper  wire 
and  are  potted  with  stycast  #1266  heat  conducting  epoxy  to  help  elim¬ 
inate  vibration  within  the  windings.  Further  design  considerations 
are  included  in  Appendix  A.  Figure  4  shows  a  close  up  photograph  of 
the  drive  coils. 

In  the  actual  installation,  the  drive  coils  form  a  part  of  a 
tuned,  series,  RLC  circuit.  The  impedence  of  this  circuit  is  given 
by: 


|Z|  =  {R2  +  (aiL  -  (wC)"1 )2}^ 


(6) 


In  this  equation  w  is  the  radial  frequency  of  the  drive  field,  L  is 
the  inductance  of  the  coils,  R  is  the  dc  resistance  of  the  coils  and 
sense  resistor  and  C  is  the  capacitance  of  the  tuning  capacitor.  In 
order  to  facilitate  tuning  and  to  allow  for  operation  over  a  range  of 
frequencies,  a  capacitance  bank  consisting  of  switch  selectable  par¬ 
allel  capacitors  has  been  constructed.  The  maximum  capacitance  avail¬ 
able  is  110  yf  which  allows  a  minimum  impedence  operation  at  frequencies 
as  low  as  200  Hz. 

The  current  through  the  coils  and  thus  the  ac  magnetic  field  is 
monitored  by  the  ac  voltage  drop  across  a  75  W,  .033  n,  "non-inductive" 
resistor.  The  inductance  of  the  sense  resistor  is  small  but  not  zero. 
However,  at  constant  frequency,  the  inductive  component  of  the  im¬ 
pedence  is  constant  and  the  potential  across  the  sense  resistor  is 
proportional  to  field  strength. 

The  sense  resistor  potential  to  field  strength  calibration  was 
calculated  by  monitering  the  induced  emf  (due  to  the  drive  field  only) 
in  the  sample  pick  up  coil.  Knowing  the  drive  field  frequency  and  ef¬ 
fective  area  of  the  coll,  the  sinosoidal  field  amplitude  can  be  related 
to  the  emf  amplitude. 


7 


The  dc  characteristics  of  the  drive  coils  are  given  in  Table 
I.  The  field  strength  was  measured  with  a  Hall  effect  gaussmeter. 

The  field  homogeniety  was  measured  to  be  better  than  1%  across  the 
surface  of  a  1"  diameter  wafer  located  at  the  central  plane  of  the 
coils. 

b.  Sample  Mounting  and  Pick  up  Coils 

The  sample  mounting  system  was  designed  for  use  with  1"  diameter 
wafers.  This  unit,  constructed  from  teflon,  has  a  1"  diameter  re¬ 
cess  in  which  the  sample  is  firmly  held.  The  holder  is  mounted  on  a 
precision  X-Y  translation  stage  which  allows  accurate  and  easy 
sample  positioning.  This  holder  can  easily  be  modified  to  hold  sam¬ 
ples  of  arbitrary  diameter. 

The  pick  up  coils  have  an  outside  diameter  of  .256"  and  were 
wound  with  #34  gauge  copper  wire  on  delrin  forms.  Both  coils  are 
mounted  on  the  drive  coil  axis  with  the  sample  coil  resting  on  the 
sample.  The  compensation  coil  is  situated  at  the  inner  edge  of  one 
of  the  drive  coils. 

Details  of  the  pickup  coil  placement  are  shown  in  Figure  5.  A 
wafer  in  the  sample  holder  mounted  between  the  drive  coils  are  also 
shown  In  Figure  5.  The  X-Y  translation  device  is  shown  in  Figure  4. 

c.  Signal  Processing 

In  order  to  null  out  the  background  signal  due  to  H^t),  a 
multi-arm  compensation  bridge  was  constructed.  This  bridge  or¬ 
iginally  designed  by  Hubert  for  use  In  M-H  loopers,  has  several 
branches  which  can  be  adjusted  to  eliminate  the  fundamental  and 
higher  harmonic  signals  from  the  pick  up  coils. 


compensation  coils 


I 

Operational  amplifiers  amplify  the  sample  and  compensation  coil 

signals  (x  10  )  before  entering  the  bridge  circuit  and  the  output  of 

2 

the  bridge  (x  10  )  before  entering  the  calculator  system. 

Signal  processing  is  accomplished  using  a  Hewlett-Packard  Model 
9825A  calculator  system.  Signal  input  is  through  a  calculator  con¬ 
trolled  HP  3437A  system  voltmeter  and  plotting  is  on  an  HP  9872A 
plotter.  The  signal  is  portioned  into  a  500  channel  array  and  signal 
averaging  is  achieved  by  sampling  each  channel  20  times.  Final 
nulling  of  the  drive  field  signal  is  achieved  by  sampling  with  a  wafer 
in  place  and  without,  and  then  subtracting  arrays.  The  array  is 
numerically  Integrated  with  respect  to  time  and  plotted  against  H^. 

An  M-H  trace  for  an  8.3  um  thick  bubble  material  is  shown  in 
Figure  6.  Traces  of  dM/dH  vs.  H  are  also  available  and  an  example 
of  this  Is  shown  in  Figure  7.  A  listing  of  the  program  used  is  in 
Appendix  b. 

d.  Calibration 

The  field  emanated  from  the  sample,  H$,  is  detected  in  the  pick 
up  coil.  Generally  in  hysteresis  loop  devices  this  field  is  pro¬ 
portional  to  some  fraction  of  the  sample  moment  m,  or  H$  =  Km.  In 
a  vibrating  sample  magnetometer,  VSM,  the  distance  from  sample  to 
detection  coils  is  large  compared  to  the  sample  dimensions.  For  this 
reason  the  constant  K  has  the  same  value  for  both  film  and  substrate. 

A  known  change  in  moment  of  the  paramagnetic  substrate  can  be  used  to 
calibrate  the  moment  axis  of  the  hysteresis  loop. 

Unfortunately,  In  the  M-H  looper  this  is  not  the  case.  The 
sample  pick  up  coll  rests  on  the  film  and  the  film  to  coil  and  sub¬ 
strate  to  coil  distances  (at  least  for  a  portion  of  the  substrate) 
are  significantly  different.  The  constants  and  Kj  are  different 
and  a  change  In  moment  In  the  substrate  (per  unit  time)  will  cause 
a  potlentlal  In  the  pick  up  coll  unlike  that  caused  by  the  same 
change  In  moment  In  the  film. 
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Fig.  6  An  M-H  trace  for  an  8.3ym  thick  bubble  film. 
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It  can  be  seen  from  the  above  discussion  that  the  gadolinium 
gallium  garnet  substrate  cannot  be  used  for  calibration  in  the 
M-H  looper.  The  use  of  a  known  or  standard  sample  is  substituted. 

The  zero  magnetization  susceptibility  of  this  sample  is  measured  using 
a  VSM. 

With  the  use  of  a  standard  (an  8.3  urn  thick  bubble  wafer)  the 
zero  magnetization  susceptibility  of  a  5.19  urn  and  a  5.6  urn  wafer 
were  measured.  As  a  check  \  was  also  measured  with  a  VSM  for  all 
samples.  These  measurements  are  given  in  table  2.  Errors  of  22%  and 
26%  resulted  in  the  two  trials. 

3.  Conclusions  and  Final  Comments 

The  high  drive  field  slew  rates  at  which  the  looper  operates 
have  certain  undesirable  side  effects.  M-H  traces  exhibit  distor¬ 
tion,  which  grows  with  slew  rate.  The  apparant  stripe-in  field, 
coercive  field  and  to  a  lesser  extent,  stripe  collapse  field  are 
strongly  influenced  by  slew  rate.  For  sample  5804  the  coercive  field 
has  been  observed  to  vary  from  3  to  40  Oe.  The  stripe  in  field 
varies  from  10  to  100  Oe. 

Delays  in  domain  nucleation  on  the  order  of  a  few  milliseconds 
and  vibration  of  the  drive  field  coils  are  believed  to  be  the  sources 
of  this  distortion.  Slight  movement  of  the  sample  and  compensation 
coils  (which  are  mounted  on  the  drive  field  coils)  result  in  an 
additional  errant,  coherent  signal  due  to  compensation  bridge  im¬ 
balance. 

At  a  drive  field  frequency  of  196.3  Hz,  x  quickly  reaches  a 
limiting  value  with  decreasing  slew  rate  for  all  samples  tested. 
However,  H  was  still  changing  (although  approaching  the  VSM  value) 
as  the  drive  field  amplitude  became  too  small  to  saturate  the  sample. 
Decreasing  the  drive  field  frequency  by  an  order  of  magnitude  would 
be  a  desirable  modification  of  this  device  and  allow  measurement  of 
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In  the  device's  present  form  noise  considerations  prohibit  the 
use  of  samples  with  thickness  less  than  about  5  pm.  Although  sucess- 
ful  loops  from  2  pm  thick  wafers  have  been  recorded,  the  extremely 
noisey  environment  at  WPAFB  does  not  allow  consistent  performance. 
More  careful  shielding  techniques,  more  sophisticated  electronics 
and  additional  stiffening  of  the  drive  coil  structure  can  overcome 
this  limitation.  Improvement  in  error  figures  should  also  be  a 
result. 

The  information  gained  in  this  research  was  presented  at  the 
24th  Annual  Conference  on  Magnetisms  and  Magnetic  Materials  at 
Cleveland,  14-18  November  1978  and  published  as  "Two  devices  for 
non-optical  characterization  of  full  wafer  magnetic  bubble  films" 

J.  Appl.  Physics  50  2182  (1979).  The  second  device  reviewed 

briefly  in  this  paper  is  the  microwave  cavity  technique  for 
measuring* bubble  collapse  and  bubble  Stripe  out  fields.  Details 
of  the  systems  are  discussed  in  the  previous  report  #AFAL-TR-78-67. 


B.  Spinwave  Excitations  in  Domain  Walls 


During  the  past  year,  the  theory  explaining  the  absorption  of 
microwave  energy  by  unsaturated  magnetic  bubble  materials  was  further 
developed.  The  work  reported  here  is  an  extension  of  a  paper  pre¬ 
sented^  and  of  the  previous  year's  annual  report^.  The  propa¬ 
gation  of  spinwaves  through  a  domain  wall  has  previously  been  exten¬ 
sively  considered^"*^.  In  this  analysis  a  single  isolated  180° 

Bloch  wall  is  considered  in  an  infinite,  lossless  medium,  with  no 
external  fields.  The  energy  density  consists  of  a  bulk  term,  limited 
to  the  uniaxial  anisotropy,  given  by 

Eb  =  -K  sin-’o  cos:>4>  .  (7) 

The  demagnetization  energy  is  ignored  initially.  The  spherical  coor¬ 
dinate  system  used  is  such  that  the  r  direction  lies  parallel  to  the 
equilibrium  magneti zation  as  shown  in  Figure  8. 

For  this  system,  the  equation  of  motion  for  the  magnetization, 

neglecting  loss  terms,  can  be  written  as 

* 

l  4  -  tb  ♦  tl  x  V'rt  (8) 

V  .It  M  • 

1  o 


where  the  first  term  on  the  right  hand  side  is  the  torque  acting  on  the 
magnetization  due  to  causes  other  than  exchange.  In  general  this  may 
include  such  terms  as  the  applied  magnetic  field,  the  crystalline 
stress  and  growth  induced  anisotropy,  and  magnetostatic  effects.  The 
second  term  is  due  to  the  exchange  interaction  between  spins.  The 
bulk  torque  term  can  be  written  as 


fb 


x  VEL 


(9) 


Here  Eb 


is  the  bulk  energy  density  excluding  exchange  effects. 


In  the  position  dependant  coordinate  system  used  here  the  total  magnet¬ 
ization  is  written  as 


M(y,t)  =  Mo  r(y)  +  my(y,t)o(y)  +  m  (y,t)*(y) 


(10) 


It  is  assumed  that  m^<<  and  m^M^.  In  this  approximation 
(M2  -  m2  -  m2  )la  s  M  . 

O  O  (j)  '  ”  o 


(11) 


For  the  coordinate  system  in  Figure  8,  a  Bloch  wall  has  4>=<j>=0.  Then  to 
first  order  in  m,  and  m., 

0  <p 

ft  x  v2ft  =  (HMe")r-  [M  me(e'  )2  +  M  m  ' '  ]e  +  [-M  mn(6‘  )2  + 

xO(|>  o  o  4>  oO 


M;>  e' '  +  M2  m' '  +  M  m..  (e‘  )2  ]<|> 

o  o  o  6 v  '  J  r  ? 


(12) 


where  the  primes  represent  differentiation  with  respect  to  y. 

The  equilibrium  configuration  is  easily  obtained  by  taking  Eqn. 

12  to  the  zeroeth  order  in  m,  and  nr  and  equating  the  right  hand  side 

(111  0  ♦ 
of  Eqn.  8  to  zero'  ' .  The  result  is 


o(n)  =  gd(n)  =  2tan  1  (en)  -  j  , 


(13) 


/  K 

where  n  is  the  reduced  variable y  and  gd(n)  is  the  Gudermanian 
function  of  n.  o(n)  is  shown  in  Figure  9.  The  following  identities 
between  o  and  n  also  exist, 


cos  o  =  sech  1 1 


or 


sin  o  =  tanh  n. 


(14) 


Considering  small  deviations  of  the  magnetization  about  the  equil- 

ibruim  position  given  by  Eq.  13,  the  equations  of  motion  for  m  and  m 

(12)  0  * 

can  be  derived'  . 

7  *  ■  eJ  *  -Jj-  [m^  cote  *  csc*9)  * 


(E°  cscecote)]  -  [(e  Y  +  mj1] 


o  '  $ 


(15a) 
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and 


—  -  ~n~  l-m  +  m  (E^  cotocsco  -  E^  cscg)]  +  rp  9 

>  ,'t  M  l'  1'0  >J>  0^  0 

For  the  energy  density  In  Eq.  7,  the  derivative  terns  are 

Eb  =  2KU  sino  cosO, 

Eb  ■  0 

* 

Ek  ■ 0 

Eoo  *  "2Ku  cos  2,1 


(16) 


e;*  =  2Ku  s1^° 

Substituting  these  into  Eq.  15  and  using  Eq.  14,  the  equations  of 
motion  reduce  to 


l  ,w 
7  ’>t 


2K, 


fT  1(1  -  2sech-  t|)m^  -  m#). 


and 


1  ^  [(1  .  2sech-’ti)m  -  m  ) , 


where  m0  represents  the  double  derivative  with  respect  to  n, 

m,  • 

0  i>rr 

These  two  equations  can  be  decoupled  by  defining  new  variables 

m,  "  m ,  •  ii\) 

0  y 


(19) 


The  decoupled  equations  are 


and 


i  >*m  _  2K„  , ,,  „  ... 

7  It  '  it  [(J  ’  2sech-  n)m+  -  m+l 


1  1(1  -  2sech-’n)m_  -  m_l 


(20a) 


(20b) 


22 


Since  both  m.  and  nr  are  real,  m*  =  m  and  Eqs.  20  are  complex 

v  $  '  " 

conjugates  of  each  other.  Assuming  a  time  dependance  of  the  form 

m+(t}  =  m±  e‘1wt  ,  (21) 


the  m+  equation  can  be  written  as 


a2ti4 
3n 


+  [n  -  1  +  2sech2n]m+  _  Q 


(22) 


where  si  =  .  Eqn.  22  is  a  Schrodinger  like  equation  with  a 

tyiSj  p 

"potential"  -2  sech  n,  which  has  a  minimum  at  n  =  0  and  goes  to  zero 
assymptotically  as  n  +  ±“.  This  is  shown  as  a  function  of  n  in  Figure 
10. 

Eqn.  22  can  be  solved  exactly  with  the  solutions  divided 

into  two  parts: 


1)  The  Bound  States:  The  bound  state  occurs  at  si  =  0  with 
the  eigenfunction  ma  =  sechn.  Since  mn  is  equivalent  to  a  change 
in  o,  one  can  write 


mQ  «  6e  =  6  sechn 
where  6  is  a  proportionality  constant.  Then 

tanhn  =  sin(e  +  6© )  =  sine  +  secose 
or 

sine  *  tanhn  -  <5ecoso 
=  tanhn  -  5  sech2n 
=  tanh(n  -  5) 

which  is  a  domain  wall  centered  at  n  =  $•  Thus  the  bound  state  corre¬ 
sponds  to  a  translation  of  the  wall.  For  the  general  potential, 

-Q  sech2n*  there  is  only  the  one  bound  state  at  si  =  0  for  Q  =  2.  If 
Q  were  to  be  increased  by  a  perturbation,  another  bound  state  would 
appear  for  n  slightly  less  than  1,  this  being  very  close  to  the  con¬ 
tinuum  states. 
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Effective  Potential 


1  and  the  eigen- 


2)  The  Continuum  States:  These  appear  for  ' 
function  can  be  written  as 


"e 


k  +  itanhn 

(1  *  k:)S 


(23) 


where  k  » ii- 1  and  the  real  and  imaginary  parts  of  Eqn.  23  are  two 
linearly  independant  solutions  of  Eqn.  22.  The  factor  of  .T+kr 
normalizes  the  solution  so  that  far  away  from  the  wall  solution  takes 
the  form 


m  =  cos  (kn  +  a) 


where  a  is  a  phase.  This  is  a  propagating  spin  wave  valid  for  both 
m.  and  m  .  The  potential  is  completely  transparent^9,10^.  However, 
the  domain  wall  distorts  the  propagating  spin  wave  as  shown  in  Figure 
11  and  develops  a  net  r.f.  magnetic  moment.  This  r.f.  moment  couples 
with  the  microwave  magnetic  field,  leading  to  absorption  of  energy. 

To  calculate  this  coupling  strength,  it  is  necessary  to  normalize  the 
solutions  to  obtain  the  amplitude  of  the  progagating  spin  wave. 

Before  normalizing,  the  inferences  of  the  demagnetization  energy  will 
be  considered. 

( g\ 

The  demagnetization  energy  is  calculated  in  the  Winter'' '  approx¬ 
imation.  The  demagnetizing  field  is  defined  by  the  following  equations 

v  •  Hd  =  =  o  (C4) 


Assuming  slow  spatial  variation  in  the  x  and  z  directions,  Eqn.  24  is 
approximately  satisfied  by 


Hd  =  -4nMy* 


(25) 


The  demagnetization  energy  then  has  the  form 
E^  *  2ntty  -  2*ffsin-'osin-> 


(26) 
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Fig.  11  Symmetric  part  of  the  spin-wave  solution  to  Eq.  (16). 


26 


Adding  this  to  Eqn.  7  leaves  all  energy  derivatives  in  Eqn.  16  un¬ 
changed  except  for  E*\  which  becomes 

<p4> 


E^  =  2Kl)sin2e  +  4iTM2sin2e 

(27) 

The  equations  of  motion 

now  have  the  form 

1  ame  ,  2IC, 

y  at  "FT 

0 

[(1  -  2sech2n)m(f  -  m^]  AttM^ 

(28) 

and 

1  3me  _ 
y  at 

(d  -  2sech2n)mQ  -  mQ] 

0 

(29) 

These  two  equations  do  not  readily  decouple.  Assuming  a  time  dependance 
e1a)t  and  drawing  on  the  experience  with  Eqn.  22,  the  following  solu¬ 
tions  are  substituted  into  Eqns.  28  and  29. 

me  =  B(k  +  itanhn)eikn 

and 

%  =  +  itanhn)e^n 

m0  and  m^  are  the  real  parts  of  Eq.  30.  This  leads  to 
y  B  »  ^  [1  +  k2  + 

and 

r c  -  ♦  *2jb 

0 

To  have  non-trivial  values  of  B  and  C 

a2  -  (1  +  k2  -  a)(i  +  k2), 

where 

» ■  At  ind  • 2  ? ' 
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(31a) 

(31b) 

(32) 


(30a) 

(30b) 


L 


a 


(33) 


Eqn.  32  can  be  solved  for  k  giving 


k-  =  \  [-2  -  a  ♦  (a2  +  4il2  )'2J 

Only  the  upper  sign  is  retained  since  when  a  =  0  (i.e.  when  demagnet- 

2 

ization  is  ignored)  k  =  s?  -  1.  Thus  the  effect  of  the  demagnet¬ 
ization  energy  is  to  change  the  value  of  k  for  a  given  frequency. 

For  typical  bubble  films,  a  -v  0.1^^.  Thus  for  q  >  1  the  change  is 
quite  small.  To  calculate  the  coupling  of  these  spin  waves  to  an 
r.f.  field  applied  in  the  xy  plane  (the  plane  of  the  film),  and 
normal  to  the  domain  wall  (hr  ^  // y),  the  effect  of  the  r.f.  field 
is  added  to  Eqns.  28  and  29.  To  first  order  in  m0  and  hy  (assuming 
h  <<  ) ,  Eqn.  29  in  unchanged  and  Eqn.  28  has  the  following  in- 

homogeneous  form. 

—  1(1  -  2sectr!n)m  -ml-  4nM  m  +  M  h 

1  "  4>  i(>  o  4,  o  y 


(34) 


A  solution  of  the  form 

m  =  ^  B.  (kccskn  -  tanhn  sinkn)  =  -  Bki*-(k,n)  (35 

0  k  K  k 

is  assumed.  The  allowed  values  of  k  are  discussed  below.  This  solu¬ 
tion  is  symmetric  about  the  origin.  The  other  form  of  the  solution 
(  k  sin  k  n  +  tanh  n  cos  k  n)  is  antisymmetric  and  will  not  give 
rise  to  a  net  moment  since  the  inhomogeneous  term  in  Eqn.  34  is 
symmetric,  m^  would  be  given  as  a  sum  in  terms  of  C^,  related  to 
Bfc  by  Eqn.  3l. 

In  oraer  to  evaluate  B^,  certain  boundary  conditions  are  re¬ 
quired  to  make  solutions  belonging  to  different  k  values  orthogonal. 
The  condition  m'  =  0  at  n  *  ±a,  where  a  »  5,  the  wall  thickness,  is 
used.  Since  the  spin  waves  are  essentially  "free"  far  away  from  the 
wall,  the  boundary  conditions  are  only  a  mathematical  convenience 
facilitating  the  calculation  of  the  k  dependence  of  the  r.f.  moment. 
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The  present  choice  insures  that  J  m^dn  equals  the  total  moment  gener¬ 
ated  at  the  wall.  The  actual  value  of  a  is  not  predicted  in  this  model 
and  was  used  as  a  fitting  parameter.  Due  to  the  presence  of  damping 
in  the  sample  the  spin  waves  excited  at  one  domain  wall  decay  before 
reaching  the  next  wall  (except  for  small  k).  Thus  the  domain  walls 
are  essentially  isolated.  This  is  achieved,  in  this  model,  by  the 
boundary  conditions  and  the  value  of  a  giving  the  best  fit  reflects 
the  amount  of  damping  present.  This  boundary  condition  puts  restric¬ 
tion  on  the  allowed  value  of  k: 


k  sin  ka  +  cos  ka  =  0. 


(36) 


With  these  boundary  conditions,  the  solutions  of  Eqn.  34  belonging  to 
various  values  of  k  are  orthogonal  (self  adjoint  differential  equa¬ 
tion  with  Neuman  boundary  conditions).  Eqn.  35  and  a  similar  relation 
for  m  ,  together  with  Eqns.  34  and  31  lead  to 


Bk  1(1  ♦  k-  +  .0(1  +  k') 


»*‘lv(k,n) 


2K 


(37) 


where  j' 

(1  +  K  + 
namely 


corresponds  to  the  applied  field  frequency.  Setting 
2 

a) ( 1  +  K  )  =  ,  and  using  the  orthogonality  of  the  * s , 


(38) 


where  both  k  and  k'  satisfy  Eqn.  36,  the  following  relation  for  is 
obtained. 


8k  " 


2y 


(39) 


V 
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A  similar  relation  can  be  written  for  C^,  using  Eqn.  31. 

The  component  of  the  moment  in  the  r.f.  field  direction  associated 
with  a  sample  of  thickness  t  and  area  is 


M 


d*w  * 


m 


hrf  dv  * 


C.iJ>(k,n)dn  •  ^s— 


-a 


2AgiMohy 


hs  ni'-ny 


[  /  v(k,n)dn]- 

-a _ 

a 

[  j  l|r’  (K,n)dr,] 


(40) 


iilkj] 
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where  is  the  stripe  period  at  zero  bias  field.  This  factor  is 

necessary  due  to  the  periodicity  of  the  domain  walls  present.  In 

the  experiment  this  moment  is  compared  with  the  r.f.  moment  of  the 

( 14) 

uniform  mode  for  which 


M 


2a  Y  M  „  ,  hy 


(41) 


The  power  absorbed  by  the  sample  i;i  each  node  is  proaortional  to  the 
cort esponding  moment.  Hence  the  ratio  of  the  power  absorptions  is 
given  by 


2  [  /  >l’(k,n)dt  ] 

-  Jldw  -  -a _  l+k: 

Mu  a  ‘  p 

P,[/  .  (k,-  )d  ] 

-  d 


(42) 


Here  has  been  set  equal  to  s:  .  This  does  not  indicate  a  resonance. 

0 

In  the  experiment  it  is  assumed  that  only  the  modes  with  the  same 
frequency  as  that  of  the  driving  r.f.  field  are  excited,  leading  to 
power  absorption  at  that  frequency;  thus  justifying  s:  =  and 
k  =  k  . 

The  integrations  in  Eq.  42  were  performed  numerically.  Similar 
calculations  for  hr<r  parallel  to  the  Z  axis  give  a  ratio  two  orders 
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of  magnitude  smaller.  Eq.  42  predicts  only  the  k  dependance  of  the 
intensity  ratio.  The  actual  magnitude  of  the  ratio  depends  on  the 
choice  of  a,  which  is  not  determined  In  this  theory.  Hence  a  is  used 
as  a  fitting  parameter. 

According  to  Eq.  40,  the  domain  wall  moment  changes  as  the  inverse 
of  the  stripe  period.  As  the  bias  field  Is  increased,  the  stripe 
pattern  changes  and  the  stripe  period,  P,  increases.  The  change  in  the 
demagnetization  field  can  be  assumed  to  cancel  the  change  in  the  bias 
field.  Thus  to  the  first  order,  it  can  be  assumed  that  the  domain  wall 
structure  remains  unchanged  and  the  only  effect  of  the  bias  field  is  to 
change  P. 

According  to  this  picture,  the  domain  wall  signal  should  change 
as  1/p  as  a  function  of  the  bias  field.  Starting  with  Eqns.  2  and  3 
of  Kooy  and  Enz.(16)the  following  relations  are  obtained. 


1 

h 


H 

4irMs 


1_  rOTT  ,  M  \ .  fnnp 

n7”  sin  ti-O  +  idH-pr 


S  -p-  s1nm"(i + 


n«l 


sinh  nms, 

~TT  1 


M  \ , rsinh  n na, 
T^)1[ - f - 1 


(43) 


(44) 


where 


i  is  the  film  thickness.  dt  and  d.<  are  the  widths  of  the  domains  with 
the  magnetization  parallel  and  antiparallel  to  the  bias  field,  respec¬ 
tively. 

Eqns.  43  and  44  are  solved  numerically  to  give  dt  +  d;  as  a  func¬ 
tion  of  II.  The  program  used  on  the  HP  9825  is  listed  in  Appendix  c 
lig.  12  shows  1/P  as  a  function  of  H,  along  with  the  integral  of  the 
actual  domain  wall  signal  for  the  sample,  Sperry  4/3/77-1.  (Actual 
signal  is  derivative  of  the  power  absorption).  The  former  was  scaled 
to  match  the  integrated  signal  at  H=0.  The  agreement  between  theory 
and  experiment  is  excellent.  Indicating  that  a  single  domain  wall 
model,  such  as  the  present  one,  is  adequate  to  explain  the  shape  of 
the  observed  low  field  signal. 

EXPERIMENTAL 

During  the  past  year,  spectrometers  operating  at  c-band  (6  GHz) 

and  L^-band  (13.5  GHz)  were  constructed.  Using  these,  data  were  taken 

on  a  number  of  samples.  FMR  data  was  taken  at  X-band  as  a  function  of 

the  angle  between  H  and  the  film  normal.  The  data  were  fitted  to 
app  /,?\ 

the  dispersion  relation  for  a  (111)  oriented  film:' 

H(.  cos ( i'-.\ )  -  241(cos  a  +  cos2,»)  -  j,  (g,(.»)  +  g.(.O)  + 

1  (45) 

|  |  { 4 n M ( cos’  .x  -  cos2.0  f  g  •  (.» )  -  gj  (,x)  )-*  +  4(^-)-'|l1 

gt(a)  =  Hulg-sin-'2o  -  ?  cos'ht  -  sin2.t(l  +  ~  sin-\0l 

7  .T 

cos4a  +  -g-  (sin2.i  -  2sin4o)l  (46) 
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Mg.  12  bolid  line:  Integral  of  the  low  field  signal  after 
subtracting  the  base  line  due  to  the  uniform  mode. 

Dotted  line:  The  Inverse  stripe  period  scaled  to 
match  the  solid  line  at  zero  field. 
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Period 


Hr  is  the  field  position  of  the  uniform  mode  and  p  is  the  angle  be¬ 
tween  H  „  and  the  film  normal.  The  sample  is  oriented  so  that  H..,. 

app  -  app 

lies  in  the  (110)  Plane.  The  angle  »,  between  M  and  the  film  normal, 
is  determined  by  the  equation 

Hsin(tf-a)  +  2t'M  sin2a  +  g,(a)  =  0  ^ 

The  data  was  fitted  to  this  relation  using  a  generalized  least  squares 
fitting  routine'  '  run  on  the  9325  computer.  A  listing  of  this  pro¬ 
gram  is  included  in  Appendix  n.  This  program  gives  the  values  of 
4iM-H  ,  and  H,  .  The  thickness  used  was  that  supplied  with  the 

U  >  K 

sample.  M  was  determined  using  the  VSM  and  A  was  obtained  through 
the  characteristic  length  determined  from  a  measurement  of  the  stripe 
period^^.  Knowing  the  material  parameters  makes  it  possible  to  com¬ 
pare  domain  wall  data  with  the  theory. 

Following  the  above  procedure  data  were  taken  on  a  number  of 
samples.  One  sample  (3-13-16)  was  studied  at  C,  X,  and  K^-band  to 
provide  experimental  data  as  a  function  of  frequency.  The  results 
ot  these  studies  were  reported  at  the  AIP  meeting  in  San  Francisco  in 
January  1978(20). 

To  obtain  more  data  on  the  frequency  dependance  implied  by  Eq.  42, 
the  domain  wall  and  uniform  mode  signals  were  studied  in  the  tempera¬ 
ture  range  180k  -  300K  at  X-band  for  the  sample  Sperry  4/3/77-1.  As 
the  temperature  is  lowered, H  increases  thus  decreasing  Thus  exper¬ 
imental  values  for  R  could  be  obtained  as  a  function  of  s.',  though  the 
actual  frequency  remained  constant.  The  temperature  dependance  of 
4 • M  was  studied  on  the  VSM  and  values  of  P  in  this  temperature  range 
were  obtained  by  using  a  cold  stage  assembly  on  the  Reichert  micro¬ 
scope.  Pi  was  found  to  be  independent  of  temperature  within  the 
accuracy  of  the  measurements. 
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In  Fig.  13,  the  reduced  ratio  is  plotted  as  a  function  of 

the  wave  vector  k.  In  Eqn.  42  the  right  hand  side  is  dependant  on 

1+k2  1 

temperature  due  to  the  factor  — - —  which  by  Eqn.  32  is  (f+'p.cO  ’ 
and  a  changes  with  temperature.  To  correct  this,  the  factor 
^ ^  is  included  in  the  L,  H,  S  and  this  reduced  ratio  is  plotted. 
Also  shown  are  the  experimental  points,  indicating  a  very  good  agree¬ 
ment  with  the  theory.  These  results  were  presented  in  the  MMM  Confer¬ 
ence  in  Cleveland  in  November^.  The  parameters  are  listed  in  Table 
3  &  4. 


TABLE  3: 

MATERIAL  PARAMETERS 


Sample:  Sperry  4/3/77-1 
Composition:  (YSmLu) 3 (FeGa) 50j z 

h(thickness) :  3.1  pm 

P  (stripe  period):  3.6  pm 
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EXPERIMENTAL  RESULTS 


Wave  Vector  k 


Fig.  13  The  reduced  Intensity  ratio  Rrt/(l+k?)  as  a  function  of 
the  wave  vector  k.  The  solid  line  shows  the  theoretical 
calculations.  The  data  points  are  shown  by  astericks. 
Data  was  taken  at  9.2  GH  .  Best  fit  was  obtained  for 
a-0.3  P0.  2 


C.  Radio  Frequency  Domain  Wall  Oscillations 


During  the  first  part  of  the  past  year,  data  was  taken  on  a  variety 
of  samples  at  a  variety  of  frequencies  using  the  RF  spectrometer  des¬ 
cribed  in  the  previous  annual  report^  Data  obtained  from  7  dif¬ 

ferent  samples  at  frequencies  from  5  MHz  to  100  MHz  show  a  low  field 
power  absorbtion  characteristic  of  the  presence  of  domains.  This 
broad  band  characteristic  indicates  that  the  effect  being  observed 
is  not  the  bound  spin  wave  state  expected  when  this  project  was  init¬ 
iated.  This  conclusion  is  further  supported  by  the  failure  of  the  cal¬ 
culations  reported  in  another  section  B  of  this  report  to  predict  bound 
states  in  this  frequency  range  (.0035  -  >:  «.  .03).  Considerable  effort 
has  been  expended  in  an  attempt  to  understand  the  nature  of  the  coupling 
of  the  in-plane  RF  magnetic  field  to  the  domains. 

In  addition  to  the  above  mentioned  work,  some  studies  of  the  operat¬ 
ing  characteristics  of  the  spectrometer  have  been  undertaken.  Of  partic¬ 
ular  interest  were  signal  strength  vs.  input  RF  power  and  signal  strength 
vs.  bridge  isolation.  As  shown  in  Figs.  14  &  15,  the  signal  strength 
varies  directly  with  RF  power  and  seems  to  be  independent  of  isolation. 
Experiments  were  also  performed  by  operating  the  spectrometer  as  a  di¬ 
rect  absorption  spectrometer,  dispensing  with  the  modulation  field  and 
lock-in  amplifier.  A  typical  spectrum  taken  in  this  mode  is  illustrated 
in  Figure  16. 

As  data  was  taken  on  a  number  of  samples,  it  became  evident  that 
it  would  be  desirable  to  be  able  to  measure  the  frequency  dependance  of 
the  domain  wall  signals.  This  is  not  practical  with  the  type  of  spec¬ 
trometer  initially  used  as  its  tuning  is  very  narrow  band  (see  Fig.  17). 
This  narrow  band  characteristic  is  a  direct  result  of  using  a  resonant 
coil  to  couple  the  rf  magnetic  field  into  the  sample,  so  the  decision 
was  made  to  modify  the  spectrometer  to  operate  in  a  non-resonant  manner 
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RF  power 


Fig.  16  Typical  direct  absorption  spectrum 


Bridge  isolation  vs.  frequency-resonator  type  spectrometer. 
Spectrometer  will  only  operate  when  isolation  is  -5dB  or  greater. 


thus  enabling  swept  frequency  measurements.  The  method  adopted  is  to 
replace  the  tuned  coil  entirely  with  a  transmission  line  having  a 
helical  center  conductor  and  terminated  in  its  characteristic  impedance. 
This  system  operates  in  the  spectrometer  in  exactly  the  same  manner  as 
the  tuned  coil,  provided  its  characteristic  impedance  was  equal  to  the 
impedance  of  the  terminator  on  the  opposite  arm  of  the  bridge.  The 
advantage  to  this  system  is  that  it  requires  no  modification  to  the 
rest  of  the  spectrometer,  except  the  addition  of  a  swept  signal  gen¬ 
erator  as  the  rf  source. 

The  design  of  the  helical  transmission  line  is  shown  in  Fig.  18 
and  is  based  on  the  characteristic  impedance  of  such  a  device,  given 
by 

z  -  £ 

7.4  t:r  (48) 

where  n  is  the  pitch  of  the  helix  in  turns/inch,  d  is  the  diameter  of  the 
helix  in  inches,  and  D  is  the  diameter  of  the  outer  conductor  in  inches^2  ^ 
Constraints  were  placed  on  the  design  by  machining  limitations.  The  val¬ 
ues  of  the  parameters  used  were  n  =  20  tums/in.,  d  =  .704  in.,  and 
D  =  .785  in.,  giving  a  designed  impedance  of  50  ft.  Copper  magnet  wire 
(28  ga)  was  used  to  form  the  helix  itself.  Due  to  assembly  difficulties, 
d  was  decreased  during  assembly  to  .684  in.,  changing  the  theoretical 
impedance  to  60  n.  After  assembly,  the  frequency  characteristics  and 
impedance  were  measured.  Using  a  pulse  technique,  the  characteristic 
impedance  of  the  helical  line  was  found  to  be  200  ft.  A  plot  of  its 
frequency  response  is  indicated  in  Figure  19.  Due  to  the  serious  mis¬ 
match  between  the  helix  and  the  rest  of  the  system,  a  resistive  impedance 
matching  network  was  designed  to  match  the  helix  to  50  ft.  This  device 
works  very  well,  and  the  response  of  the  system  with  the  helix  on  the 
bridge  is  shown  in  Figure  20. 
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Fig.  19  Frequency  response  of  helical  transmission  line 


transmi ssi 


To  date,  no  data  has  been  taken  with  this  device  and  work  is 
continuing  toward  assembly  and  use  of  a  spectrometer  using  this  de¬ 
vice. 


2:  Theory 

( r><> ) 

Drawing  on  work  done  by  Slonczewski .  the  generalized  equations 
of  motion  for  a  domain  wall  are 


AY 


(q  -  *.vi) 


^q 


(49) 


Here  q  is  the  wall  position,  s  is  the  angle  of  the  magnetization  with 
respect  to  the  x  axis,  a  is  the  Gilbert  damping  coefficient  for  the 
material,  A  is  the  wall  width  parameter  (a  y-  )  *') .  and  v*  is  the  en- 

Nu 

orgy  density  of  the  wall.  The  coordinate  system  used  is  shown  in  Mg. 

21.  The  functional  derivative  is  defined  by 

>Nq 

Ao  _  >K'  A.' 

Aq  ~  ,iq  *  avq  (SO) 

Including  terms  to  account  for  exchange,  uniaxial  anisotropy,  demag¬ 
netization  and  ground  state  energy,  the  variation  of  the  energy  density 
becomes 

i  2A 

Ao  =  A |  —  (vq)'  ♦  2Aa ( vy ) ;  Asin'O  -  2ij)|  + 

(51 ) 

"M.\(Hxsinv  -  HyCOSY ) Ay  -  2MHzAq 

where  Hx,  H^,  and  H2  are  the  cartesian  components  of  the  applied  mag¬ 
netic  field  and  v  »  x  ♦  7  . 
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Fig.  21  Fixed  coordinate  system  xyz 
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Performing  the  functional  derivatives  as  indicated  in  Eqn.  (49) 
and  rearranging  the  resulting  equations  gives 


iflq  *  2irMYSin  2(y  -  |^)  +  -^-(Hxsinf  -  Hycosv)'^p  vy  +  ay 
t  -  YHZ  -  2nMYA  ^(si n2(v  -  ff)  +  ^  ^  '  £  $ 


(52) 


To  a  first  approximation,  vq  =  vy  =  0  at  the  center  of  a  thin  film. 

Also,  if  the  x  axis  lies  in  the  plane  of  the  wall,  ^  =  0.  These  two 

a  X 

approximations  simplify  Eqn.  (52)  so  that  it  becomes 

A-1q  =  2ttM  y  sin2y  +  —  (H  sins'  -  H  cosy)  +  ay 

c  x  y 

(53) 

If  there  is  no  bias  field  (H  *  0)  and  H  =  H  »  0  at  t  3  0,  the  initial 

l  x  y 

conditions  are  q (0)  *  y(0)  =  0  and  the  equation  for  S'  can  be  solved  as  a 

function  of  q  .... 

(54) 

*(t)  -  -Jq(t) 


Thus,  Eqn.  (53)  becomes  a  single,  first  order  differential  equation 

<  ■  ’'TO  **  t  cos(i  ,)  sihUll  <55> 

Here  it  has  been  assumed  that  Hx  =  0  and  Hy  -  Hrf  sin(wt). 

Using  material  parameters  measured  in  this  laboratory  on  sample 
3-18-16  (see  table  5),  numerical  solutions  were  obtained  for  Eqn.  (55). 
A  typical  solution  is  shown  in  Fig.  22. 

Based  on  the  numerical  solutions,  analytic  solutions  to  equation 
53  were  sought.  If  Hz  -  0,  i  -  ^  q  implying  that  y  *  ~  q.  From  the 
numerical  solutions  It  Is  seen  that  j  q  is  small  enough  that  slnU'  *  v 
and  cosy  •  1.  Thus,  Eqn.  (53)  becomes 
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TABLE  5 


Material  parameters 

4  M  -  259  G 
Hu  -  1164  Oe 
A  -  3.7  x  10"2  erg/cm 
*  1.76  x  107  He/Oe 
h  -  4. 1pm 
1  =  .50pm 
a  =  3.2  x  10-2 
PQ  -  8.33pm 
H  -  95  Oe 


of  sample  3-18-16 

saturation  magnetization 
uniaxial  anisotropy  field 
exchange  constant 
gyromagnetic  ratio 
film  thickness 
characteristic  length 
Gilbert  damping  coefficient 
zero  field  stripe  period 
cubic  anisotropy  field 
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(56) 


( 1+a2  )q  =  -4irMYaq  -  —  (H^q  +  AHy) 


Solutions  to  Eqn.  (56)  can  be  found  for  the  special  cases  =  0  and 
Hy  =  0.  For  the  case  H*  =  0  and  Hy  =  H^sinwt  the  solution  is 
dU)  =  - nA-Hcf - sin(wt  +  <(>) 

2»[(5^)2  +  (1*^)]'= 

tan*  = 

4nMya 

For  the  case  H  =  o  and  H  =  H  f$inwt  the  solution  is 
y  xiT 


q(t)  =  q(0)e 


~4itMyoi  iryHrfaCOSu>t 

e  . 2uj(  l+a2  ) 


(57) 


(58) 


Eqn.  (57)  represents  an  oscillitory  motion  of  the  domain  wall,  with  a 
phase  lag  between  the  wall  motion  and  the  driving  field.  The  ampli¬ 
tude  of  the  oscillations  is  a  complicated  function  of  frequency  and 
varies  linearly  with  the  magnitude  of  the  drive  field.  Eqn.  (58),  on 
the  other  hand,  represents  an  exponential  decay  of  wall  displacement 
modulated  by  a  complicated  oscillitory  factor  at  the  frequency  of  the 
drive  field.  For  typical  materials,  the  decay  constant  is  on  the  or¬ 
der  of  108  sec-1,  so  that  compared  to  the  steady  state  behavior  of 
Eqn.  (57)  this  solution  is  only  a  transient  and  need  not  be  consider¬ 
ed  in  most  cases. 


Consider  now  the  case  where  H 
time  dependence  si  nut.  Then  Eqn. 


=  0  and  H 
*  y 
(56)  becomes 


=  0  but  both  have  the 


where 


q  +  (A  +  Bsinwt)q  =  Csinwt 


A  - 

4nMTa 

A  - 

l+a2 

B  = 

2 

C  = 

-it  HVA 

2(1+ a2 

g 

l+a2 


(59) 


(60) 


Assume  a  solution  of  the  form 


q(t)  *  u(t)qQ(t) , 

where  qQ(t)  satisfies 

qQ  +  (A  +  Bsinwt)qQ=  0. 

It  can  easily  be  shown  that  the  solution  to  Eqn.  (62)  is 

Q 

(At  +  -  cos  <ot) 

q,(t)  ■  q  (0)e 


(61) 


(62) 


(63) 


Using  this  fact,  a  differential  equation  for  u(t)  is  found  to  be 

“  ■  '-At  e  6/"  m 


whose  solution  is 

u  ( t)  -  — 1 e‘At  e  B/u>  COSult  sinmt  dt  (65) 

v’o'  ' 

Now,  consider  the  case  Hx  s  -  0,  H  »  Hsinwt.  In  this  case, 
Eqn.  (53)  becomes 


A"1  q  *  2nM y  sin2y  +  e»Y 

f  *  y  H  Sinut  -  7  q 
A 


(66) 


The  two  parts  of  Eqn.  (66)  can  be  decoupled,  and  making  a  small  angle 
aprroximation  for  v,  the  resulting  equation  is 

(l+a;)Y  +  4nM>aH'  =  YH  sinu>t 


The  solution  to  Eqn.  (67)  is 


f(t)  -  V"3t  +  la^T^T 


sin  (wt  +  <i>) 


(67) 

(68) 


where 


a  = 


b  * 


4  n  MY 

I+S7 

YH 

T+* 


(69) 


*»♦  ■  1  -  ’^> 

and  y  is  a  constant  determined  by  the  value  of  v  at  t  =  0.  In- 
serting  this  solution  into  Eqn.  (66)  gives  a  differential  equation 
for  q.  This  equation  includes  terns  proportional  to  e-at.  Since 
only  steady  state  solutions  can  give  rise  to  a  long  duration  power 
absorbtion,  these  terms  can  be  ignored.  The  equation  for  q  then 
becomes 

q  -  (^yA)((a-j  sinU  ♦  4>)  ♦  (a COs{lot  +  {?0) 


The  solution  to  Eqn.  (70)  is 

q(t)  =  C  sin(u't  +  *  +  )  +  D 

where 


(4irMY^b)->(bu)g):>A 


a-  +  u> 


(71) 

(72) 


and 


tan<j>'  =  "InJliA. 


UKX 
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i. 


i 


and  the  value  of  D  is  determined  by  the  wall  displacement  when 
sinwt  *  0. 


In  order  to  calculate  the  power  absorbed  by  the  wall  motion,  it 
is  necessary  to  know  the  force  acting  on  the  wall.  According  to 


Slonczewski 


A22) 


,  the  wall  momentum  is 
.  .  2M  ... 


and  assuming  that  H,  =  0,  the  force  acting  on  the  wall  is 

F  =  p  =  —  f 

Y  YA 

so  that  the  work  done  in  moving  the  wall  is 

.  _  -2Ma  ij. 
dw  _  _____  qdq 

The  power  absorbed  is  then 

p  -  dw  _  -2Ma  jj,  . . 

P  -  dt  "  ~  dt 

Using  the  solution  given  by  Eqn.  (57),  this  is 

- -  cos2(cot  +  *)  -  sin2 (wt  ♦  *) 

4[(*^r  4  ( 1+a2 )'  J 


(73) 

(74) 

(75) 

(76) 

(77) 


A  plot  of  the  amplitude  of  Eqn.  (77),  using  the  material  parameters  of 
sample  3-18-16,  is  shown  in  Fig.  23.  A  plot  of  this  same  relation¬ 
ship,  using  material  parameters  of  sample  TI330,  along  with  data 
taken  on  TI330,  is  shown  in  Fig.  24.  Agreement  between  experiment 
and  theory  seems  to  be  reasonable. 

Note,  too,  that  since  Hrf  a/P^n“  ,  where  P^n  is  the  power  input 
to  the  coil,  P<*Pjn*  in  good  agreement  with  the  data  of  Fig.  . 

Thus,  this  theory  seems  to  provide  an  adequate  explanation  of  the 
data,  and  in  particular  provides  a  mechanism  whereby  an  inplane  rf 
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Fig.  24  Power  absorption  vs 
Theory  curve  is  lea 
points. 
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magnetic  field  can  drive  domain  wall  motion.  More  experimental  work 
needs  to  be  done  to  evaluate  the  qualitative  and  quanitative  predic¬ 
tions  of  this  theory,  and  the  modified  spectrometer  described  in  the 
previous  section  will  be  capable  of  obtaining  the  necessary  data 
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SECTION  II 


Testing  Program 

A.  Overlap  Pattern  Verification 

Perhaps  the  largest  class  of  visual  pattern  recognition  problems 
involve  the  identification  of  objects  whose  presentation  is  either 
different  in  time  or  whose  context  is  variable.  For  example  the 
recognition  of  aircraft  (23)  involves  the  identification  of  an  ob¬ 
ject  whose  orientation  and  size  is  constantly  changing.  Face 
identification  (24,25)  has  similar  problems  plus  variations  caused 
by  changes  in  facial  expression,  makeup,  hair  style,  etc.  Many  more 
examples  could  be  cited  all  of  which  require  object  identification  in 
a  changing  environment  (26).  There  appear  to  be  few  problems  of  real 
interest  which  involve  patterns  fixed  in  space  and  time.  Certainly 
little  appears  in  the  open  literature  concerning  such  problems.  The 
advent  of  bubble  memories  (27)  however  has  created  interest  in  just 
such  a  problem. 

The  simplicity  and  low  number  of  masking  operations  required  in 
the  construction  of  a  bubble  memory,  make  possible  the  development 

9  1  ? 

of  extremely  large  memories  of  the  order  of  10  to  12  bits  (28). 
Because  of  the  immense  size  involved,  exhaustive  testing  of  each 

9 

memory  device  becomes  excessive.  Consider,  for  example,  a  10  bit 
memory  running  at  500  k  bits/sec.  Such  a  memory  requires  2000  sec.  or 
3.3  hrs.  to  access  every  bit  once.  Thus  a  complete  memory  test, 
including  tests  for  pattern  errors,  could  take  several  days.  Ob¬ 
viously  some  gross  pretesting  is  necessary  to  reduce  the  number  of 
devices  which  must  be  exhaustively  tested.  Clearly  if  a  fault  free 
magnetic  substrate  (garnet,  for  example)  and  a  fault  free  permalloy 
overlay  and  aluminized  conduction  paths  exist,  then  the  probability 
that  a  memory  device  works  Is  extremely  high.  This  sorting  on  the 
basis  of  these  two  factors  should  reduce  the  overall  testing  time  of 
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the  memory  elements.  Further,  in  order  to  increase  yield  rates,  man¬ 
ufacturers  currently  allow  some  overlay  faults  to  exist  if  these 
faults  can  be  isolated  and  removed  from  service  without  affecting 
the  rest  of  memory  (29).  Thus  any  overlay  faults  detected  should 
be  marked  so  that  specific  tests  can  be  performed  to  determine  their 
effect  on  the  rest  of  the  device.  The  purpose  of  this  paper,  there¬ 
fore,  is  to  describe  the  experimental,  optical  test  procedure  which 
has  been  shown  feasible  in  checking  overlay  patterns  for  faults  in 
four  micron  bubble,  bubble  memories. 

Basically  the  problem  involves  detecting  differences  between  two 
fixed  geometry  patterns.  One  of  these  patterns,  the  fault  free  one, 
is  taken  as  the  authority  and  the  other  is  compared  to  it  in  some 
appropriate  fashion.  Such  a  procedure  is  generally  termed  template 
matching  (30).  The  major  difficulty  with  any  pattern  recognition 
problem  is  the  selection  of  the  features  to  be  used  in  making  the 
comparisons.  It  is  tempting  in  problems  where  the  geometry  is  ab¬ 
solutely  fixed  to  use  a  pattern  description  language  (31,32)  to  de¬ 
fine  a  correct  pattern  and  then  compare  the  derived  description  of 
the  sample  to  that  of  the  authority.  Unfortunately,  deriving  a 
description  of  a  pattern  is  generally  a  time  consuming  proposition. 

It  would  seem,  therefore  that  statistical  recognition  techniques 
(33,34),  which  generally  require  less  time  for  feature  extraction 
than  do  syntactic  techniques  (35)  would  be  the  procedure  of  choice. 

The  selection  of  features  to  be  used  in  a  geometrical  or  statis¬ 
tical  approach  to  pattern  recognition  is  perhaps  the  single  most 
important  factor  in  the  success  of  recognition.  Unfortunately,  there 
is  little  in  the  way  of  a  general  theory  to  help  in  the  selection  of 
these  features.  However,  for  the  problem  being  considered  here  three 
factors  are  important: 
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1)  the  features  should  be  rapidly  computable 

2)  they  should  be  capable  of  isolating  the  location 
of  a  fault  in  the  field  of  view  and  not  merely 
identify  that  a  fault  exists,  and 

3)  only  a  small  number  of  features  should  be  required. 

The  first  requirement  is  due  to  the  very  large  number  of  views 
necessary  to  cover  a  large  scale  bubble  memory  chip  with  adequate 
resolution.  More  will  be  said  about  this  later.  The  second  require¬ 
ment  is  necessary  so  that  faults  which  affect  only  portions  of  a  mem¬ 
ory  chip  can  be  isolated  and  removed  from  service  in  some  way.  The 
last  requirement  is  necessary  so  that  computer  memory  space  is  kept  to 
some  reasonable  value. 

Autocorrelation  coefficients  (36),  Fourier  transform  coefficients 
(37)  and  moment  invariants  (38)  have  all  been  used  with  varying  de¬ 
grees  of  success  in  optical  pattern  recognition  problems.  Certainly 
criteria  1  and  3  above  are  satisfied  for  these  features.  Unfortunate¬ 
ly,  all  of  these  features  tend  to  represent  the  image  in  a  rather 
"global"  way,  that  is,  little  or  no  information  is  obtained  about 
localcharacteristics.  Histograms,  however,  appear  to  meet  all  three 
requirements  and  therefore  were  the  features  selected  for  this  prob¬ 
lem. 

The  images  used  are  binary,  black-white  images  which  are  mapped 
onto  a  128  x  128  bit  binary  array.  Two  sets  of  histograms  are  ob¬ 
tained  from  this  array;  horizontal  and  vertical.  The  ith  element 
of  the  horizontal  histogram  represents  the  number  of  white  cells  on 
row  i  of  the  image  array.  Similarly  for  the  vertical  histogram  com¬ 
ponents.  If  it  is  assumed  that  only  a  single  "blob-like"  fault 
occurs  within  the  field  of  view  then  a  comparison  of  the  horizontal 
and  vertical  histograms  with  those  of  the  authority  should  localize 
the  fault  in  the  image  plane. 
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EXPERIMENTAL  SYSTEM 

Generally  large  bubble  memory  devices  are  integrated  on  large 
(perhaps  a  few  square  inches)  substrates.  It  is  obvious  then,  that  a 
very  large  number  of  views  will  be  required  to  cover  the  device  with 

b 

adequate  resolution.  For  example  a  10  bit  memory  with  16  micron 
elements  on  a  one  square  inch  substrate  would  require  roughly  10,000 
separate  nonoverlapping  views  to  cover  the  entire  surface.*  Position¬ 
ing  the  device  under  a  microscope  with  the  accuracy  necessary  is,  of 
course,  an  extremely  difficult  problem.  However,  systems  such  as  the 
Leitz  Step  and  Repeat  camera  (39),  which  is  capable  of  positioning 
to  within  ±.2  microns  might  be  capable  of  handling  the  problem.  More 
will  be  said  about  this  in  section  B  below. 

Assuming  that  the  positioning  problem  is  solvable,  the  recognition 
task  becomes  one  of  acquiring  an  image,  preprocessing  to  clean  and 
smooth  it,  extract  the  histogram  features,  and  finally  determine  if 
faults  occur.  The  images  used  in  this  work  were  obtained  by  a  Zeis 
model  microscope  from  an  experimental  bubble  memory  device  of  approx¬ 
imately  24,000  bits.  Figure  25  shows  one  such  photograph.  These 
photographs  were  next  blown  up  to  an  8x10  size  and  digitized  by  a  TV 
sytem  connected  to  a  DEC  PDP-9  computer  (41,42).  The  resulting  bi¬ 
nary  images  were  then  smoothed  in  an  128  x  128  binary  image  array 
processor  (43)  also  connected  to  this  computer.  Figure  26  shows 
the  resulting  control  image  used  for  the  experiments  described 
in  what  follows. 

Two  experiments  were  carried  out.  The  purpose  of  the  first  was 
to  determine  the  efficacy  of  the  histogram  features  for  identification 
and  isolation  of  overlay  pattern  faults.  The  second  experiment  was 
intended  to  determine  the  effects  of  positioning  error  on  fault  de¬ 
tection.  In  the  training  phase  of  experiment  one,  the  control  image 

♦This’  assumes  an  128  x  128  image  and  that  two  pixels  per  4  micron 
line  width  are  required  for  adequate  resolution. 
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Fig.  25  Original  Overlay  Pattern. 


was  input  to  the  system  10  successive  times  without  motion  and  the 
horizontal  and  vertical  128  number  histograms  were  computed  at  each 
trial.  The  average  value  and  standard  deviation  were  then  computed 
for  each  of  the  256  histogram  components.  The  resulting  average  hist¬ 
ograms  with  their  supplementary  standard  deviation  information,  was 
then  used  as  the  authority  for  the  recognition  phase  of  the  experiment. 
This  averaging  process  was  carried  out  to  reduce  the  effects  of  camera 
chain  noise  on  the  results. 

In  the  recognition  phase  of  the  first  experiment  overlay  pattern 
faults  were  simulated  by  either  adding  to  the  pattern  or  occluding 
portions  of  the  pattern.  This  was  done  by  placing  irregularly  shaped 
white  and  black  pieces  of  paper  on  the  input  photograph.  Two  sizes 
were  used  for  these  faults;  a  large  error,  roughly  the  size  of  the  top 
of  one  of  the  T's,  and  a  small  error,  about  the  width  of  a  bar  in  the 
image.  The  modified  image  was  then  input  to  the  system  and  the  hist¬ 
ograms  computed.  For  each  component  of  the  histogram  a  binary  class¬ 
ification  number  was  computed  as 


Ci  = 


h.-H. 

1 

0  otherwise 


(78) 


where  h^  is  the  ith  component  of  the  sample  histogram,  is  the  mean 
of  the  authority  histogram  component  i,  SDi  is  the  standard  deviation 
of  the  authority  histogram  component  and  T  is  an  arbitrarily  set 
threshold.  In  order  to  avoid  problems  of  quantitization  noise  which 
occurs  when  a  scan  line  is  on  the  edge  of  image  element,  a  second 
classification  number  D..  was  introduced  as 
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D1  ■ 


1  if  C.=Ci+1-l  or  Ci_1=Ci=l 
0  otherwise 


(79) 


i  *  ,2,. . . ,255 

where  the  first  128  components  correspond  to  the  horizontal  histogram 
and  the  last  128  correspond  to  the  vertical.  Figures  27,  28,  29,  and 
30  show  the  results  of  this  experiment.  In  these  figures  the  Di  is 
plotted  along  the  edge  of  the  picture  and  used  T=6  standard  deviations. 
The  important  thing  to  observe  is  that,  indeed,  the  faults  in  the 
overlay  patterns  were  detected  and  furthermore  were  located  by  the 
positioning  of  the  values. 

In  order  to  determine  the  effects  of  positioning  error  on  the 
recognition  accuracy  a  second  experiment  was  carried  out.  In  the 
training  phase  of  this  experiment  the  control  image  was  allowed  to 
vary  by  ±1/2  pattern  bar  width  in  both  the  horizontal  and  vertical 
directions.  The  authority  histograms  were  then  computed  on  the  basis 
of  two  samples  taken  at  each  of  the  eight  possible  off  center  positions 
plus  the  centered  image.  These  18  samples  were  used  to  compute  the 
control  histogram  component  averages  and  standard  deviations. 


Equations  78  and  79  were  again  used  in  the  recognition  phase. 

The  additive  and  subtractive  faults  were  simulated  as  before  and  are 
shown  in  Figures  31  through  38.  The  plotted  in  Figures  31  through 
34  used  a  threshold,  T,  of  6  as  in  the  first  experiment  and  those  of 
Figures  35  through  38  used  a  threshold  of  2.  The  additive  faults 
shown  in  Figures  31  and  32  are,  indeed,  marked  by  the  D...  However, 
the  fault  in  Figure  32  is  detected  only  in  the  vertical  histogram 
and  not  the  horizontal.  Figure  33  shows  a  large  subtractive  fault 
in  the  horizontal,  but  not  the  vertical,  histogram.  The  subtractive 
fault  of  Figure  34  is  undetected  at  threshold  level  6. 

Figures  35  through  38,  however,  which  represent  the  same  faults 
but  use  a  threshold  of  2,  do  show  all  faults  as  being  identified. 
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fig.  33  Large  Subtractive  fault,  16 
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However,  some  false  alarms  start  to  appear.  The  obvious  Triplication 
is  that  there  is  a  coupling  between  positional  variation  of  the  de¬ 
vice,  the  optimum  threshold  for  detecting  the  faults,  and  the  rate  of 
false  alarms. 

CONCLUSIONS 

It  is  obviously  impossible  to  draw  any  major  conclusions  from  the 
small  sample  size  given  here.  However,  it  is  equally  as  clear  that  the 
histograms  do  appear  to  mark  the  faults,  both  additive  and  subtractive, 
at  least  for  the  samples  taken.  Furthermore,  it  would  seem  that  posi¬ 
tioning  variations  of  up  to  1/2  bar  .width  are  tolerable.  Thus  the  use 
of  histogram  features  for  detection  of  faults  in  fixed  geometry  prob¬ 
lems  would  appear  to  be  feasible. 

B.  Positioning  Problem 

As  stated  above,  the  general  approach  being  taken  in  the  verifica¬ 
tion  of  the  overlay  pattern  is  to  scan  the  overlay  in  several  steps  each 
time  comparing  the  acutal  image  with  a  corresponding  perfect  image. 
Errors  will  show  up  as  differences  in  the  comparision.  Since  the  over¬ 
lay  pattern  is  quite  large,  it  will,  in  general,  be  necessary  to  move 
the  pattern  relative  to  the  optical  axis. 

Because  of  this  some  mechanical  positioning  mechanism  is  required 
to  move  the  overlay  array  relative  to  the  sensor  with  at  least  two 
mutually  perpendicular  degrees  of  freedom.  (Assumes  angular  orienta¬ 
tion  already  established).  The  requirements  and  constraints  applicable 
to  each  of  these  degrees  of  freedom  are  nearly  identical  so  it  will 
suffice  to  discuss  only  one. 


It  is  generally  true  that  mechanical  positioning  mechanisms  are 


frought  with  numerous  inherent  nonlinear  characteristics  which  result 
in  position  errors  which  are  rarely  repeatable  from  one  trial  to  the 
next.  Dead-zone  due  to  gear  train  backlash  in  the  presence  of  fric¬ 
tion  probably  constitutes  the  major  source  of  error.  The  extent  to 
which  nonlinear  effects  can  be  masked  is  a  function  of  the  type  of 
servomechanism  used. 

There  are  two  basic  types  of  electrically  actuated  positioners. 
First,  the  open  loop  or  open  cycle  type.  Typically,  these  consist 
of  a  step  motor  driving  a  load  through  a  low  backlash  gear  train.  In 
such  a  drive  the  accuracy  of  positioning  of  the  load  is  determined 
by  the  dynamic  and  static  characteristics  of  the  step  motor  and  by 
the  backlash  and  frictional  properties  of  the  gear  train  as  well  as 
the  friction  and  inertia  of  the  load. 

A  typical  step  motor  indexes  200  steps  per  revolution  with  non- 
accumulative  angular  position  error.  However,  the  static  error  (a 
quantity  not  usually  specified  by  manufactuerers)  at  a  given  step 
position  must  be  considered.  The  torque  available  to  the  load  is 
approximately  a  linear  function  of  angular  position  deviation  about 
the  desired  position,  is 

T  ~  itla  360° 

T  =  K( 0  -  200®  ^ 

where  e  is  the  actual  angular  position  in  degrees 

n  Is  the  increment  number  for  a  200  step  per  revolution  motor 
K  is  the  steady  state  torque  constant 

The  magnitude  of  K  is  a  function  of  the  size  of  the  motor  and  the 

magnitude  of  the  excitation  current. 

The  dynamic  characteristics  of  step  motors  must  also  be  consid¬ 
ered.  An  unloaded  step  motor  is  governed  by  an  underdamped  second 
order  characteristic  equation,  i.e. ,  its  angular  shaft  position  leaves 
a  steady-state,  approaches  and  overshoots  the  new  desired  position, 
then  oscillates  with  an  exponentially  damped  amplitude  about  the  new 


steady  state  position  and  finally  stops.  A  mechanical  load  in  the 
shaft  alters  the  overshoot,  ringing  frequency  and  rate  of  exponential 
damping.  Nearly  always  an  unrealistically  large  amount  of  external 
viscous  damping  is  required  to  eliminate  overshooting.  (The  excep¬ 
tion  is  that  one  can  decrease  overshoot  for  a  given  load  by  decreas¬ 
ing  motor  current.  The  sacrifice  is  usually  intolerable  since  K  and 
steady  state  position  error  increase  with  decreasing  excitation  cur¬ 
rent.)  In  some  instances,  it  is  possible  to  use  some  type  of  Geneva 
mechanism  to  mask  the  effects  of  overshoot  and  ringing. 

Thus  a  very  careful  mechanical  design  must  be  incorporated  in 
step  motor-driven  open-loop  designs  to  maintain  error  at  an  accept¬ 
able  level. 

The  second  basic  type  of  electrically  actuated  posi toner  is  a 
closed-loop  type  in  which  a  motor  (not  a  step  motor)  overcomes  fric¬ 
tion  and  moves  through  the  backlash  deadzone  until  the  load  reaches 
the  desired  position.  At  such  a  time  the  motor  exerts  whatever  torque 
is  necessary  to  maintain  the  load  at  the  desired  position  within  a 
position  error  of  ±  e.  In  this  case  however,  given  a  type  1  servo 
controller,  the  magnitude  of  e  is  primarily  a  function  of  the  resolu¬ 
tion  limits  of  the  transducer  sensing  the  position  of  the  load.  In 
the  extreme  case  it  is  conceivable  that  one  could  use  an  optical  de¬ 
fraction  grating  positon-sensing  transducer,  e.g.,  the  Leite  Step  Re¬ 
peat  Camera  in  which  a  ±  .2pm  position  accuracy  is  specified.  (The 
width  of  permalloy  overlay  bars  in  the  sample  available  is  approximate¬ 
ly  3um). 

Obviously  the  closed-loop  approach  involves  a  higher  degree  of 
complexity  and  expense  however,  position  accuracy  and  repeatability  far 
exceed  those  of  an  open-loop  configuration. 

The  following  consists  of  a  discussion  of  positioning  accuracy 
requirements  for  two  broad  classes  of  patterns:  a)  those  patterns 
which  are  repeated  in  both  the  x  and  y  directions,  and  b)  non-repeated 


patterns  consisting  of  isolated  features  or  long  narrow  strips  of 
features.  It  is  also  to  be  expected  that  positioning  accuracy  will 
be  strongly  influenced  by  the  limitations  of  the  pattern  recognition 
technique  used. 

First,  it  should  be  recognized  that  a  closed-loop  positioner  with 
an  accuracy  of  ±  .2utn  opens  options  not  available  with  lower  accuracy 
positioners.  For  example,  no  overlapping  of  scanning  would  be 
necessary,  software  search  modes  to  locate  a  non-periodic  feature 
would  not  be  required,  the  position  of  the  boundary  of  a  repeated  pat¬ 
tern  would  be  precisely  known.  An  attempt  will  be  made  to  specify 
the  extent  to  which  lesser  accuracies  would  be  feasible. 

The  type  of  pattern  imposing  the  least  stringent  positioning 
accuracy  demands  is  one  which  is  repeated  or  periodic  in  both  X  and 
Y.  Such  a  pattern  could  be  examined  in  successive  narrow  overlapping 
strips,  each  string  being  divided  into  n  overlapping  subsections.  The 
data  collected  from  each  subsection  would  be  compared  with  idealized 
data  from  a  model  pattern.  Obviously  if  there  are  p  basic  X-cycles 
and  r  basic  Y-cycles  in  the  mode,  the  portion  of  sensed  data  actually 
used  must  contain  exactly  p  and  r  basic  X  and  Y  cycles  respectively. 
There  are  two  important  conclusions  from  the  foregoing.  First,  sensor 
deflection  manitude  and  rate  must  be  extremely  closely  controlled 
since  it  is  the  scan  velocity  which  bridges  the  gap  between  actual 
length  and  the  apparent  length  of  a  scanned  object.  It  is  conceivable 
that  a  closed-loop  scan  control  using  an  edge  groticule  sensed  by  the 
scanner  as  feedback  would  be  required. 

The  second  conclusion  is  that  since  only  a  portion  of  the  data 
obtained  from  scanning  a  subsection  of  the  overlay  pattern  can  be  used, 
the  overlap  of  the  scanned  subsections  must  be  determined  from  the 
dimensions  of  the  unused  portions  of  the  subsection.  Specifically, 
the  distance  between  center  of  scanned  subsections,  D,  must  be 
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D  <  P  .  5x  -  2e  -  6x  =  (p-1)  5x  -  2e 

where  6x  is  the  X  dimension  of  a  basic  X-cycle 

P  is  the  number  of  basic  X  cycles  in  the  model 
±e  is  the  maximum  position  error. 

From  this  inequality  one  can  conclude  that  for  a  repeated  pattern 
nearly  any  position  error  would  be  acceptable  provided  the  increased 
processing  time  resulting  from  greater  overlap  could  be  tolerated. 

As  long  as  processed  areas  did  not  include  the  boundary  of  a 
periodic  array,  the  scheme  described  appears  to  be  viable.  Special 
problems  arise  when  boundaries  are  inevitably  reached.  A  means  would 
have  to  be  developed  for  distinguishing  between  sensed  data  discrep¬ 
ancies  from  the  ideal  due  to  faulty  patterns  and  that  due  to  the  scan¬ 
ning  of  a  legitimately  incomplete  subsection. 

Non-repeated  portions  of  a  pattern  may  pose  more  stringent  pos¬ 
itioning  accuracy  requirements  than  the  repeated  type.  If  relatively 
low  positioning  accuracy,  i.e.,  error  greater  than  one  basic  cycle 
were  considered,  some  type  of  search  mode  would  be  required  to  locate 
an  isolated  narrow  line  of  overlay  patterns.  The  searching  and  sub¬ 
sequent  processing  of  such  a  pattern  would  have  to  be  capable  of  dis¬ 
criminating  between  erroneously  deposited  blobs  and  legitimate  pat¬ 
terns. 

Another  type  of  problem  will  exist  in  the  data  collection  process 
regardless  of  the  pattern  configuration.  A  given  scan  line  will  in¬ 
evitably  run  along  the  black-white  or  1-0  boundary  of  a  section  of 
overlay.  The  resulting  raw  data  will  be  a  somewhat  random  sequence 
of  l's  and  0's.  The  question  which  must  be  addressed  is,  can  such 
noise  be  successfully  handles  without  excessive  compromise  or  must 
provisions  be  made  to  make  minute  scan  position  corrections  when  such 
a  situation  is  detected? 
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If  the  need  to  examine  the  basic  cycles  in  the  vicinity  of  the 
edge  of  a  continuous  array  did  not  exist,  and  also  if  no  non-periodic 
patterns  required  processing,  postioning  accuracies  well  within  the 
capability  of  open  loop  step  motor  positioners  could  be  used.  How¬ 
ever,  since  such  patterns  must  be  included  the  question  remains:  can 
sufficiently  sophisticated  pattern  recognition  techniques  be  incor¬ 
porated  to  compensate  for  the  rather  gross  position  uncertainties  of 
an  open  loop  positioner,  or  must  the  positioner  be  made  to  be  more 
sophisticated  to  eliminate  position  ambiguity? 

Another  requirement  completely  independent  of  the  foregoing  is 
that  scan  velocity  must  be  sufficiently  closely  controlled  to  avoid 
time  varying  magnification  errors. 

Unfortunately,  do  to  early  termination  of  this  contract,  little 
experimental  work  was  carried  out  to  determine  requirements  on  posi¬ 
tioning  accuracy.  Some  preliminary  results  have  been  presented  in 
part  A  above  and  indicate  the  feasibility  of  the  approach. 

C.  Optical  Waveguide  Work 

Once  the  correctness  of  the  overlay  pattern  is  verified,  dynamic 
testing  of  the  chip  is  required.  In  order  to  speed  this  testing  up  it 
is  necessary  that  multiple  detection  paths  be  present.  As  proposed 
here  these  paths  are  optical  in  nature  and  generally  require  that 
polarized  light  be  coupled  into  the  substrate  region  and  then  coupled 
out  of  this  region  after  having  the  polarization  shifted  by  Faraday 
Rotation  due  to  any  bubbles  in  the  light  path. 

The  experiments  performed  in  the  ear’y  months  of  the  contract 
were  directed  toward  coupling  light  into  thin  film  waveguides  by 
using  gratings.  This  method  was  used  to  aid  in  the  physical  contact 
between  the  magnetic  bubble  memory  chip  to  be  tested  and  the  integrat¬ 
ed  optic  test  circuit.  Two  different  approaches  were  considered  while 


using  gratings.  Figure  39  shows  the  light  coupling  into  the  film  from 
the  air,  while  Figure  40  shows  the  second  method  of  coupling  the  light 
from  the  substrate  side.  These  two  methods  were  evaluated  by  first, 
spinning  on  KTP  positive  photoresist  at  300  RPM  and  contact  printing 
a  grating  of  .1  mil  lines  spaced  300  mil  apart  using  standard  photo- 
liography  procedures  on  a  YI6  film  grown  on  a  GG6  substrate.  Next 
the  sample  was  placed  in  an  experimental  apparatus  so  that  a  HeNe 
laser  could  be  directed  toward  the  front  as  shown  in  Figure  39  and 
the  back  as  shown  in  Figure  40.  The  intensity  of  the  input  was 
measured  to  be  .78  mwatt  while  the  maximized  output  was  measured  to  be 
0.34  mwatt  for  the  front  and  .063  for  the  back.  The  output  was  max¬ 
imized  in  each  case  by  varying  the  angle  of  incident.  Since  coupling 
from  the  front  has  no  mechanical  advantage  and  the  coupling  efficiency 
is  less  than  that  of  coupling  from  the  back  side,  it  leads  to  the  con¬ 
clusion  that  coupling  from  the  back  is  a  better  approach.  However 
10*!.  coupling  efficiency  is  not  good  enough  since  the  lost  energy  is 
scattered  throughout  the  device  and  surrounding  space.  Another  prob¬ 
lem  arose  when  the  magnetic  moment  parallel  to  the  surface  was  re¬ 
fracted  in  a  localized  area  by  applying  an  external  magnet.  When 
measuring  the  output  the  percentage  change  in  the  light  level  was 
too  small  to  detect.  From  these  results  it  was  concluded  that  a  dif¬ 
ferent  method  of  coupling  light  into  the  waveguide  was  needed. 

By  using  rutile  prisms  which  are  highly  birefingent  a  TM  wave 
being  decoupled  from  a  thin  film  waveguide  will  be  separated  by  20oll' 
from  the  TE  wave  emerging  from  the  same  field.  Since  the  Faraday  ro¬ 
tation  effect  of  the  magnetic  substrate,  ^FgOj;?,  rotates  the  polar¬ 
ization  by  280  /cm,  propagation  of  the  light  through  lOnmeters  will 
give  a  6  rotation  of  .28°.  Using  a  rutile  prism  will  lend  well  in 
detecting  such  a  small  ao. 

The  experimental  appratus,  as  shown  in  Figure  41,  will  consist 
of  Y3Fe5®l2  substrate.  The  rutile  prisms  will  be  used  to  couple  light 
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into  and  out  of  the  YIG  epilayer.  This  device  will  be  in  the  rotating 
magnetic  field  when  in  close  proximity  of  a  magnetic  bubble  memory 
chip.  However  the  sampling  time  of  the  detection  should  be  when  the 
magnetic  field  is  in  such  a  position  as  to  minimize  the  TE  mode  while 
no  bubble  is  in  the  detectors  presence.  This  will  result  in  maximum 
sensitivity  when  a  bubble  is  detected.  With  this  setup  and  using  a 
1  miliwatt  1.152  um  He-Ne  laser  a  delta  intensity  of  approximately 
10  uwatts  will  be  detected,  which  is  much  better  than  previous  results. 

As  discussed  before,  prisms  were  not  used  at  first  due  to  mechan¬ 
ical  problems.  With  the  added  detecting  capability  it  is  hoped  prism 
coupling  will  prove  to  be  more  advantageous  to  use. 


F/ 8  9/2 
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APPENDIX  A 


Helmholtz  Coils 

The  drive  colls  used  In  the  M-H  looper  are  of  the  Helmholtz  de¬ 
sign.  Helmholtz  colls  consist  of  two  colls  arranged  such  that  the 
distance  between  the  coll  planes  Is  equal  to  the  mean  radius  of  the 
colls.  In  such  an  arrangement,  the  field  In  the  central  portion  of 
the  median  plane  Is  fairly  homogeneous. 

A1 

The  components  of  the  field  between  the  colls  are  given  by  : 

Hz  -  ^  [0. 44959  (1  -  g^r)  +  0.0035967  (jlb^36a2_)  j£  (3  COs<’o-l) 
-  0.06474  (35  cos40  -  30  cos^O  +  3)] 


and 

Hx  "  W  s  1  nocoso [ - •  0035967 ( 3-6a‘)  +  .06474  {jj-  (28  cos-'o 

-12)1.  (Al) 

The  coordinate  axis,  r,  o,  and  a  and  b  are  defined  in  Figure  Al.  I  is 
In  amps,  n  is  the  number  of  turns  in  each  coil  and  R,  a,  b  and  r  are 
in  cm.  The  equation  assumes  that  the  current  is  distributed  uniformly 
over  the  crossection  a,  b. 

In  the  median  plane  between  the  coil,  o  =  n/2.  In  this  case; 

Hx-° 

and 

HZ  *  1 r  1,44959  (1  “  SSlH  *  •  0035967 ^-)  -  3(.06474)j£) 

(A2) 

Figures  A2  and  A3  are  a  listing  and  typical  output  of  a  program  (for 
a  H-P  9820  calculator)  to  plot  the  magnitude  of  the  field  In  the 
median  plane.  In  figure  A3,  the  coordinate  r  represents  the  distance 
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Figure  A-2.  Program  for  HP  9820  calculator  used  to  calculate  the 
magnetic  field  of  the  M-H  looper  drive  colls.  At 
the  right  Is  a  typical  printout. 


Figure  A-3.  A  plot  of  the  magnetic  field  In  the  Median  plane. 


from  the  central  point  In  the  median  plane.  Hmax  is  the  value  of  the 
field  at  the  center.  The  program  calculated  that  for  these  dimensions, 
H/Hmax  *  when  r  =  2,83  cm*  Theoretically  the  field  should  be 
homogeneous  to  1%  over  a  2  Inch  diameter  circle.  The  parameters  used 
in  this  calculation  were  those  used  for  the  final  coils. 


APPENDIX  B 


Signal  Processing  Program  for  the  M-H  Looper 

The  computer  program  (for  use  with  a  HP  9825A  calculator) 
used  for  signal  processing  on  the  M-H  looper  is  listed  below. 

0:  rad;fxd  3 

1:  dim  A*  (3J,  B$  (3],  X  [500],  Y  [500],  Z  [500],  T  [500],  Y  [500] 
2:  time  1000 
3:  on  err  "PLOTTER" 

4:  rds  (705)  -  T 
5:  on  err  "DMV" 

6:  rds  (725)  -  Z 
7:  time  0 
8:  wrt  705,  "IN" 

9:  wrt  725,  "N1SR2T2F1" 

10:  dsp  "COMMAND  ?";  beep;stp 

11:  "PLOTTER":  dsp  "PLOTTER  DOWN";  stp;  gto  3 

12:  "DVM" :  dsp  "DYM  DOWN";  stp;  gto  5 

13:  "A":  0-*-A  +  M->-F  +  S--P  +  G;ent  "AVER  AGING?:,  A$ 

14:  ent  "FREQUENCY?",  F 
15:  ent  "SINGLE  SWEEP?",  B$ 

16:  if  B$  =  "yes";  1  -  M 
17:  ina  X,  Y,  Z,  Y,  T 
18:  for  I  *  1  to  500 
19:  I/(500*F)  •+■  T [ I ] 

20:  next  I 

21:  for  I  -  1  to  500 

22:  fmt  "D",  f.7,  "S";  wrt  725 

23:  red  725,  Y[I];  dsp  I;  next  I 

24:  if  A$  »  "no';  gto  32 

25:  for  0  »  1  to  19 

26:  dsp  J 
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27:  for  I  *  1  to  500 

28:  fmt  "D",  f.7,  "S";  wrt  725 

29:  0  •>  G;  red  725,  G;  Y[I];  +  G  -  Yfl];  next  I 

30:  next  J 

31:  smpy  .05*Y  -  Y 

32:  if  M  =  1;  9to  36 

33:  ara  Y  •*  Z;  1  -*  N 

34:  beep;  dsp  "INSERT  SAMPLE";  stp 

35:  gto  21 

36:  beep;  ent  "CURRENT  AMPLITUDE  (Amps)?",  A;  A-24,  59 

37:  for  I  =  1  to  500 

38:  Acos  (2"FT[Il)  *  XlH 

39:  next  I 

40:  gto  10 

41:  "0":  ina  V;  0  *•  S;  0  -*•  P;  Tor  I  =  1  to  500 
42:  Y[ll  >  S  *-S;  2[IJ  +  P-*-P;  S  -  P  ^Vll) 

43:  next  I 

44:  min  (V[*l)  *  C;  max  (VI*])  -  E 

45:  scl  -A,  A,  C,  E;  pen#  3 

46:  for  I  *  1  to  500 

47:  pit  XII]*  VIII;  next  I 

48:  pen 

49:  gto  10 

50:  "G":  scl  -A,  A,  C,  E;  pen#  4 

51:  xax  0,  A/5,  -A,  A,  5 
52:  yax  0,  -C/10,  C,  E,  5;  pen 
53:  gto  10 

54:  "C":  ina  V;  ara  Y  -  Z  -*■  V 

55:  for  I  *  11  to  239 

56:  V{l]/ab$(sin(2wFT[I]))  -*•  V[I];  next  I 

57:  for  I  «  261  to  489 

58:  Vll]/abs(sin(2nFTtI]))  -*•  VII];  next  I 


A;  ina  X 
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59 

for 

I  - 

490  to  500 

60 

0  - 

V[I] 

;  next  I 

61 

for 

I  » 

1  to  10 

62 

0  - 

V[I] 

;  next  I 

63 

for 

I  • 

240  to  260 

64 

0  - 

V[I1 

;  next  I 

65 

min(V(*] )  -  C;  max(V{* 

1) 

66 

scl 

-A. 

A,  C,  E;  pen# 

3 

67 

for 

I  - 

1  to  500 

63 

pit 

X(I1 

,  V[ I ] ;  next 

I 

69 

pen 

70 

gto 

10  *12610 

APPENDIX  C 


Program  to  calculate  stripe  period  as  a  function  of  the  bias  field. 


0  :  rods  sf  g  14s snp  "#  of  points".  Us  dim  ALU,  21s  fit 

1  t  snp  "Ku".H.  “4nM".G.“P0".P.  “Thick".  T 

2  i  8nH/G-Hs 1+G/H-rls frl-r7s 0-K 

3  «  P/Tr7-P»  P-Xs  0-Qs gsb  "EVAL" 

4  i  -Y-Ks  prt  "1/h*", Kj spo 

5  «  . 9/U-Gs P-Ai  1-J 

6  i  "STP3"s 30A-Bs  P/30-Ds  ls-9-Es 20-M 

7  ,  prt  "M/Ms-".  Qt gsb  "ROOTS" 

8  i  0-Os  for  N*1  to  100 

9  i  o  i  n  (Nit  ( 1 +0)  )  -r2»  sxp  (Nn/X)  -r3 

1 0i  r3- 1 / r3-r4»  r3+ 1 / r3-r5t  r4+r7r5-r6 

lit  Y+(Xr2r4/NNr6-r8)-Ys nsxt  N 

12s  "SER"t  sin  (Nn  (1+0)  )  -r2s  sxp  (Nir/X)  -r3 

1 3s  r3~ 1 /r3-r 4s  r3+ 1 /r3-r5s  r4+r7r5-r6 

14s  Y+(Xr2r4/NNr6-r8)-Y 

15s  if  abs (r8) >abs (Y) ls~9s  N+l-Ns  gto  “SER" 

16s  (2r7Y/ imi+Q) *G-A LJ.  11s  Xr7-ACJ.  21s  prt  "H/4irMs=".  ALJ.  11s  spo  3s  J+l-J 
17,  if  (Q+.  9/U-Q)  <■•  9s  X-As  gto  "STP3" 

18s  dsp  "Plottsr  on  ?"s stp 

19s  sol  0.  A  CU.  11 . 0.  1/P S  pan#  Is  xax  0.  A  CU.  1)  /Us  yax  0. .  05s  1-Os  pan#  2 
20s  pit  A  CO,  11. 1/A  CO.  21s  pans  if  (Q+l-Q) <°Us jmp  0 
21s  snp  "Stors  inFils#".  Ws ref  W,  AL*1 
22s  and 

23s  "R00TS"s  ofg  2s  if  A>Bssfg  2s  dsp  "ILLEGAL  BOUNDS^  stp 
24s  A-Xsgsb  "EVAL" 

25s  Y-F 

26s  if  (A+O-A)  >Bs  dsp  "No  Root  Found"sstp 
27s  A-Xsgsb  "EVAL" 

28s  if  (FY-Z)  >0s  gto  -3 
29s  if  Z<8s  gto  +3 
30s  if  F«0s A-D-Xs  F-Y 
31s  A+D-As ls-12-Ss gto  +9 
32s  A-O-Ls  A-Rs  0-C 
33s  (L+R) /2-Xs  gsb  "EVAL" 

34s  if  (C+l-C)>Ms  gto  +5 
35s  if  abs(Y)<Es0-Zsgto +2 
36s  if  (FY-Z»0s  X-Ls  gto  -3 
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37*  if  Z*0*  R-L-S*  gto  +3 
38*  X-Ri  gto  -5 
39*  R-L-S 

40*  (R+D/2-X*  prt  "P/h«",  Xr7* apo  iprt  "ACCURATE  TO".  Stapc  i  rot 

41*  "EVAL"*0-Y 

42*  for  N*1  to  20 

43*  ain  (N*  (1+0) /2)  ~*r2t  axp  (Nir/X)  "*r3 

44*  r3**l/ r3“*r4*  r3+l/r3-*r5;  r4+r7r5-*r6 

45*  Y+(r2r2(Nar7/r6-Xr4)/NNNr6-^8)-Y*naxt  N 

46*  "SUM"*ain(Nir(l+0)/2)-r2iaxp(Nir/X)-r3 

47*  r3-l/r3-*r4j  r3+l/r3**r5j  r4+r7r5**r6 

48*  Y+  (r2r2  (Nnr7/r6*-Xr4)  /NNNr6-r8)  -Y 

49*  if  aba (r8) >aba (Y)  le-9*  N+l“*Ni  gto  "SUM" 

50*  K+2rlXY/*nnr*Y*rot 
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APPENDIX  D 


Curve  flttinR  proRram  for  FMR  data. 


0  ,  »nt  "Nuafcor  of  doto  poiro".  N.  "Nuobor  of  poraootoro*.  M 
1  ,  dip  ACW.BCW.CCH3.OtM3.FtW.GCH.M3.HCW.KCN3 
2.  dip  ICW.JCM.W.LCN3.SCW.HCW.XCW.YCW 

3  i  dim  AtC33 

4  i  oot  "Uoing  unoortaintioo?".  At 

5  i  if  At-"yoo"i  ofg  1 

6  .  prt  "Xi"."Yi"iif  flgliprt  "Sigpa  i“ 

7  i  opo 

0  t  for  1*1  to  N 
9  ,  ont  "Xi"»XCI3,  "Yi".  YCI3 
10i  if  flgliont  "Sigpa  i\SCI3 

11.  prt  XC13. YCI3.  if  flgliprt  SCI3 

12.  opo 

13.  noxt  I 

14.  opo 

15s  prt  "Initial  paranotoro"i  opo 

16s  for  1*1  to  Mi  ont  "Parapotor".  ACI3|  prt  ACI3ino*t  I 
17.  opo  i oppy  .  BIA^D 
10.  .  001-Fi  oil  'Corfu' 

19s  C-Ai  prt  A 

20.  "Loopotart"i  ol 1  'Curfit 

21.  if  <C~A) /A<*.  01i  gto  "Oono" 

22i  C-Ai  prt  A 

23.  gto  "Loopatart" 

24.  "Dono"* 

25.  prt  Ci  opo 

26.  prt  "Parasiotor",  "Sigs»o"i  opo 

27|  for  1-1  to  Mi  prt  ACI3.  CCI3i  opo  inoxt  I 

2ft  prt  "fit  valuoo"i  opo 

29.  for  1*1  to  Niprt  FCIlinoxt  I 

30i  opo  3 

31s  ond 

32i  "Curfit"i 

3ft  N-M-Piif  P<*ft  ft»Ci  rot 

34.  ino  Vi  l.H.  G»Li  1 

3Si  if  flgliara  S*S-*Ki  ara  L/1MI 

3ft  for  1*1  to  N 


ii 
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37*  oil  'Fdoriv' 

38t  for  J-l  to  M 

39k  HCJl^m*(YCI3-  ,Funotn,(I»*UJ]41U] 

40i  for  K-l  to  J 

41*  GCJ.K3MIU3MU3MCK3-GU.K3 

42*  noxt  K 

43*  noxt  J 

44*  noxt  I 

45*  for  J-l  to  Mi  for  K"1  to  J|GU.K3-GCK.  J3inoxt  Kinoxt  J 

46*  for  IM  to  N 

47*  'Functn' (I)-FCI3*noxt  I 

40*  'Fohiocj'-D 

49*  "71"*  for  >1  to  M 

50k  for  K-l  to  M 

51*  GCJ,K3/f (GCJ.  J3GCK.K3WU.K3 

52*  noxt  K 

53*  l+F-JCJ.  J3 

54*  noxt  J 

55*  inv  J-J 

56*  for  J-l  to  H 

57,  ACJ3-BCJ3 

58,  for  K-l  to  M 

59*  BLJ3+H  CK3  *  J  [  J,  K3  /{  (G  CJ.  J3  G  CK.  K3  >  -B  CJ3 

60*  noxt  K 

61*  noxt  J 

62*  for  1-1  to  N 

63,  ‘Funotn’CD-FCHinoxt  I 

64*  ’Fohioq'-C 

65*  if  D<C*10F-F,gto  "71" 

66,  for  J-l  to  M 

67*  BCJ3-ACJ3 

60*  f (JCJ, J3/GCJ, J33-CCJ3 

69*  noxt  J 

70,  .  1F-F 

71*  rot 

72*  "Fdoriv", 

73*  for  J-l  to  M 


74*  A  [J]  -Xi  D  CJ3  -Zi  X+Z-A  CJ3 
75*  'Functn'CD-Y 
76*  X-Z-A  LJ] 

77»  CY“  *  Fonotn '  (I) )  /ZL*\  [J] 

78*  X-ACJ] 

79*  noxt  J 

80i  rot 
81*  "Fohiocj"* 

82*  0-C 

83i  for  I“1  to  N 

84*  owm*(Ym-Fm)t2-c 

85i  noxt  I 

86  i  rot  C/P 

87  t  ond 

88  t  "Funotn"* 

89  *  0-G-Bj  *f g  2j  XCplJ-0 

90  i  "TWO"*  YCI]oin(0-B)*.5AC2]  <oin<2B)**rl)**E 

91  .  ACl](l/24*rl(l+7(ooo(2B)-r2))-f2/6*oin(B)t2(1^2r2))^E-€ 

92  .  if  ob*(D  <10-91^0  "OUT" 

93  «  if  EG<0igto  "SIX- 

94  *  B-LiE-G»B+O/10-B-Rigto  "TWO"*  if  B>90»  d*p  "Q"»  0-B*  gto  "OUT- 

95  «  "SIX"*  <L*R)/2-B 

96*  YCI]#in(0-B)+.5AC2]  (oin(2B)-*rl)**E 

97  «  Am(l/24*rl(l+7(eo#(2B)-r2))-f2/6*oin(B)t2(l*2r2))*E-E 

98  *  if  abo (D <lo”9» gto  "OUT- 

99  *  if  EG<0*B-Rigto  "SIX"iif  R-L<.001,  <R*L)/2-B»gto  "OUT- 
100*  B-L»E-G*gto  "SIX"*  if  R-L<.001i<R*L)/2-ftgto  "OUT- 
101*  -0UT-«*in(B)t2-p2»0-B-KCpl] 

102*  oin<2B>-p3 
103*  I“p2-*p4 
104*  p4-p2-*p5 
105*  coo(4B)-p6 

106*  A  Cl]  (3/16+35p6/48+5p5/12+5f2p3p2/3)-*p7 
107*  A  C13  (G-7p6+4p5)/16+f2p3(5p2/3-l))'*p6 
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